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Abstract 

The beyond mean field approximation combined with the Gogny interaction is used 
to analyze the spherical to prolate deformed shape transition in the Neodymium 
isotopic chain. The vibrator as well as the rotor limits are nicely reproduced while 
the transitional region is only qualitatively described probably due to the lack of 
triaxial correlations in the calculations. Our results do not support the interpretation 
of 150 Nd as a critical point nucleus and question the interpretation of shape changes 
as nuclear shape phase transitions. 
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Phase transitions have always attracted the attention of scientist in all research 
fields. In nuclear physics, in particular, shape transitions and the super-fluid 
normal-fluid transition have been extensively studied since a long time at finite 
temperature (1) in many different theoretical approaches as well as experimen- 
tally. Also at zero temperature the phenomena of shape coexistence and shape 
transition have been a fascinating working field and many papers have been 
devoted to it (2; 3). Most of these works have concentrated on the shape tran- 
sitions from spherical to prolate (oblate) or triaxial deformed shapes. In these 
studies there have been used several approaches from the mean field approach 
to the Generator Coordinate Method (GCM) with schematic interactions (4) 
and effective ones (5) as well as 'exact' calculations with small configuration 
spaces and realistic interactions (6). However, the overwhelming contribution 
to this topic has been made within the framework of the interacting boson 
models (IBM) (3; 7; 8; 9; 10; 11). In this algebraic approach, spherical, prolate 
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and 7-soft nuclei are associated to the dynamical symmetries of the model, 
U(5), SU(3), 0(6), respectively. Also the general aspects of the shape transi- 
tions have been described within the classical limit of the IBM which is the 
Bohr (collective) Hamiltonian (3). 

However, the concept of shape phase transitions is, at least, controversial. In 
finite systems one does not expect real phase transitions because the large 
quantum fluctuations may spoil the sharp behavior of the corresponding pa- 
rameters found in a macroscopic system leading to a blurred phase transition. 
In addition, the difficulties associated with the identification (non-uniqueness) 
of the order parameter as well as the lack of a continuous control parameter 
(Z, N are integer numbers) have contributed to create a certain reluctance 
of part of the nuclear community with respect to this issue. With this back- 
ground, some recent developments by Iachello (12; 13) came out as a surprise, 
since he went even further by proposing the existence of critical point nuclei 
in the nuclear shape phase transitions. Iachello developed analytical solutions 
to the Bohr collective Hamiltonian with potentials which correspond, within 
the framework of the model, to the critical point of the first order shape phase 
transitions between spherical and prolate nuclei (the so-called X(5) symme- 
try), as well as to the critical point of the second order phase transition from 
spherical to 7-soft nuclei (identified as the E(5) symmetry). Furthermore, these 
new solutions were used to assign a critical character to some nuclei which ex- 
perimentally show similar excitation spectra as the ones predicted by these 
theoretical approaches. Hence, the isotopes 152 Sm, 150 Nd were associated to 
the X(5) symmetry (14; 15) while the nucleus 134 Ba would correspond to the 
E(5) symmetry (16). 

Although the IBM and the Bohr collective Hamiltonian are useful representa- 
tions for describing shape phase transitions, it is important to mention that 
both models oversimplify the underlying nuclear many-body problem. Hence, 
the interacting boson models have many free parameters which are specifi- 
cally adjusted to a particular region of nuclei and, therefore, their predictive 
power is limited. Concerning the collective Hamiltonian, one misses a com- 
pelling justification of the potentials used in such approaches. For instance a 
derivation based on microscopic calculations with reasonable nucleon-nucleon 
interactions has not been presented yet. Furthermore, these potentials ignore 
nuclear degrees of freedom which might be significant in a transitional region 
and focus only on the quadrupole deformation. 

The aim of the present work is to analyze these shape transitions from a 
microscopic viewpoint using state of the art Beyond Mean Field Approaches 
(BMFA). In particular we will study the spherical to prolate deformed tran- 
sition in the Neodymium isotopes {Z = 60), with the intention, after com- 
parison with the experiment and the X(5) predictions, to shed some light on 
this delicate problem of the nuclear phase transitions. There have been some 
self-consistent mean field approaches with effective phenomenological interac- 
tions which have also addressed the shape phase transitions in the A ~ 150 
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region (17; 18; 19). Recently, Niksic and collaborators have reported the first 
BMFA study of the proposed critical behavior of the 150 Nd isotope (20). Their 
study is based on the Relativistic Mean Field approach and concentrates on 
the nuclei 148 - 15 °- 152 Nd. Our BMFA is based on the Hartree-Fock-Bogoliubov 
(HFB) approach and in the calculations we use the Gogny force. In this Letter 
we extend the calculations of ref. (20) to describe the vibrational and rota- 
tional limits of the shape transitions, i.e. the nuclei 144 - 146 - 148 - 15 °- 152 - 154 Ncl ; 
and we furthermore analyze the relevance of triaxial shapes along the tran- 
sition region. In our BMFA (21) we first generate the collective subspace of 
HFB wave functions using <f-constrained particle number projection before the 
variation, i.e., we minimize the energy 



F N.z,*_ (* N ' Z mH\* N > Z (q)) m 

with \<& N ' z (q)) = P N P 2 *\(f)(q)) , where P N ( Z ^ j s the projector onto neutron (pro- 
ton) number, \<fi(q)) are HFB-type wave functions and q = qi, ...,qjM is a set 
of intrinsic constraints that allows the definition of potential energy surfaces 
along the most relevant degrees of freedom. The configuration mixing calcula- 
tion is performed within the generator coordinate method (GCM) framework 
taking linear combinations of the particle number and angular momentum 
projected wave functions obtained in the first step 



Then, the variational principle applied to the weights f N ' z > J > a (q) gives the 
generalized eigenvalue problem (Hill- Wheeler equation): 



J (H N ' Z ' J {q, q*) - E N > z > J '°N N > z > J (q, <f )) f N ^( ( f)d ( f = 0, (3) 

with 7-t N i z > J and J\f N,z,J the norm and Hamiltonian overlaps, respectively, see 
(22) for further details. We use in our calculations the finite range density de- 
pendent Gogny force with the D1S parametrization (23), which is well known 
for its successful predictions. 

Let us now turn to the IBM description. The Bohr Hamiltonian allows the 
description of the quantum many-body system in terms of a one-body Hamil- 
tonian with a potential energy V which depends on the collective variables 
7) (3; 24), i.e. H co u = T(/3, ■y) + V(P, 7). The collective model can be solved 
analytically under certain symmetry restrictions and the comparison of these 
solutions with the experimental data is used to identify the intrinsic shape 
of a particular nucleus. The potential energy surfaces (PES) V(f3, 7) which 
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Fig. 1. (a)-(f) Potential energy V(/3, 7) (Eq. 4) given by the collective model to 
describe the U(5)-SU(3) transition for different values of the control parameter £. 
Contour lines are plotted each 0.05 arbitrary units, the energy is set to zero at the 
minimum of each surface. 

describe the spherical to prolate deformed transition (U(5)— ►SU(3)) can be 
obtained in the framework of the interacting boson model using the method 
of intrinsic states (13): 



N(3 2 

T+J 2 



N(N-l) 
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A(3 2 + 2V2/3 3 cos3 7 + ^/? 4 



(4) 



where N is the number of valence pairs (we take N = 10 as in Ref. (13)) 
and £, (0 < £ < 1), is the continuous control parameter which determines the 
final structure of the potential and is used to simulate the phase transition. 
In Fig. 1 we show this potential, in the (J3, 7) plane, for the most significant 
values of (. For ( = 0, see panel (a), we find a potential very stiff in the 
(3, 7 variables which correspond to a spherical nucleus. For higher ( values 
(panels (b), (c) and (d)) we find a large softening in the (3 deformation and 
a smaller one in the 7 direction. For larger ( values (panels (e) and (f)) the 
nucleus becomes prolate deformed and rather stiff in the 7 direction. One gets 
more insight looking at the axially symmetric trajectory (see Fig. 2). Here, we 
observe the evolution from a spherical potential (£ = 0), whose minimum is 
found at (3 = 0, to potentials with well-deformed minima (£ = 0.028,0.030). 
For intermediate ( values one finds a set of potential energy curves which 
are practically degenerated along the prolate axis in the interval (3 G [0,0.4]. 
These curves show two minima, an spherical and a prolate deformed one, and 
can be found only in a very reduced range of values for the control parameter 
(0.0251 < C < 0.0255). In particular, for ( = 0.02547 the spherical and the 
prolate deformed minima are degenerated and this condition defines precisely 
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Fig. 2. Potential energy (Eq. 4) without normalization along the axial trajectory 
7 = 0°,60°. 

the critical point of a first order phase transition where the order parameter is 
the deformation (3. We denote V crit . (/?, 7) the corresponding critical potential. 
The narrowness of the mentioned ( interval strongly contrasts with the integer 
character of the control parameter in real nuclei. 

To compare with the experimental data one should consider the dynamical de- 
grees of freedom, i.e., substitute V crit XP,l) in the Bohr Hamiltonian and solve 
the Schrodinger equation H co u& = E<& to find energies and wave functions 
to calculate transition probabilities. Comparison of the theoretical predictions 
with the experimental data should help to identify the hypothetical critical 
point nuclei. This program has been conducted by Iachello (13) in a very el- 
egant and clever way approximating VcritXfl, 7) in a convenient manner. He 
assumed that the potential was separable in the variables (/?, 7) and approx- 
imated the (3 dependence by a square well in the interval (3 G [0,/?m] and 
the 7 one by a harmonic oscillator. With these approximations one obtains 
analytical solutions to the Bohr Hamiltonian which describe the critical point 
of a phase transition from a spherical to an axially symmetric (prolate) shape 
and are associated to the so-called X(5) symmetry (13). 

Let us now analyze the PES of the isotopes 144_154 Nd calculated with the 
Gogny force to see how much they resemble those of Fig. 1. Since the po- 
tential of Eq. 4 preserves the number of particles and the Bohr Hamiltonian 
recovers the rotational invariance, it is obvious that at the microscopic level 
the potential V((3, 7) of Eq. 4 is given by E N ' z (q), see Eq. 1, with q= ((3,j). 
Furthermore collective hamiltonians similar to the Bohr Hamiltonian can be 
obtained from the Hill- Wheeler equation considering the Gaussian overlap 
approximation and also from adiabatic time-dependent Hartree-Fock approx- 
imations (24; 4). Therefore, the diagonalization of the Bohr Hamiltonian is 
comparable to the solution of a Hill- Wheeler equation. In Fig. 3 we show the 
PES calculated using Eq. 1, here we observe the shape transition from a prac- 
tically spherical nucleus ( 144 Nd) to well- prolate deformed nuclei ( 152_154 Nd). 
In these limits we appreciate a similarity with the potential surfaces given 
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Fig. 3. (a)-(f ) Particle number projected potential energy surfaces in the (/?, 7) plane 
(Eq. 1) calculated with the Gogny D1S force for gQ 4_154 Nd. Contour lines represent 
an step of 1 MeV and the energy is set to zero at the minimum of each surface. 

by the geometrical model (panel 1(b) in the first case and panel 1(f) in the 
deformed one). 

However, we also notice several differences between the surfaces given by the 
microscopic and the geometrical model. Thus the microscopic ones present 
more substructure due to the underlying shell structure which is ignored in 
the collective model. In addition, there is no continuous control parameter 
in the isotopic chain. Nevertheless, the most important difference is the role 
played by the triaxial degree of freedom in the calculations with the Gogny 
force for the transitional isotopes 146_150 Nd. In these cases, we notice that the 
surfaces soften not only in the axial prolate direction, as in Fig. (1), but also in 
the 7 direction. Finally, the 152_154 Nd nuclei show symmetric minima isolated 
from the rest of the (/?, 7) plane with high barriers. 

In Fig. 4 we show the potential energy curves along the axial deformation. 
Here, we observe that prolate and oblate deformed minima are obtained for all 
nuclei with a maximum at /3 = 0, in contrast to the result shown in fig. (2). The 
deformation corresponding to those minima is displaced to higher values with 
increasing mass number. Furthermore, both the energy difference between the 
prolate and oblate minima, and also the value of the barrier between them, 
are larger with larger number of neutrons in the system. Nevertheless, it is 
important to note that in the 144_150 Nd isotopes, the oblate minima are saddle 
points in the (/?, 7) plane, while for the 152_154 Nd nuclei they are real minima. 

The study of the particle number projected PES only gives a qualitative 
picture of the shape transition. To obtain spectroscopic information and to 
include configuration mixing, we perform GCM calculations along the axial 
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Fig. 4. Particle number projected potential energy surfaces along the axial trajectory 
7 = 0°,60° calculated with the Gogny D1S force for ^ 4_154 Nd. The zero of the 
energy is chosen at (3 = 0. 




Fig. 5. Collective wave functions along the axial quadrupole deformation (/?) for 
J = 0, J = 2 and J = 4, calculated for 144_154 Nd isotopes. \G((3)\ 2 represents the 
probability to find the deformation (3 in the collective wave functions, see also (22) 

quadrupole degree of freedom with particle number and angular momentum 
projected wave functions. As we have already mentioned in the previous para- 
graph, the triaxial effects may play an important role in the transitional nuclei. 
However, configuration mixing calculations that include the full triaxial an- 
gular momentum projection are not feasible with the current computational 
capabilities. Therefore, a significant deviation of the theoretical prediction 
from the experimental data only for the soft triaxial nuclei could point out to 
the relevance of this degree of freedom. 

In Fig. 5 we analyze the effects of the angular momentum projection and the 
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Fig. 6. (a) Ratio between the energies of the first 4 + , 2 + excited states, (b) re- 
duced quadrupole transition probability between the first 2 + , ground states for the 
144_154 Nd isotopes. Dots and circles correspond to experimental — extracted from 
Refs. (25; 15) — and theoretical values, respectively. 

GCM on the ground and first excite states (J 77 = + ,2 + ,4 + ). The collective 
wave functions (Eq. 2) of the ground states of the 144_146 Nd isotopes exhibit a 
noticeable configuration mixing of oblate and prolate shapes, the 146 Nd nucleus 
being more deformed than the 144 Nd one. The prolate-oblate mixing tends to 
disappear with increasing value of the angular momentum where clearly only 
prolate mixing is obtained. On the other hand, concerning the 148_154 Nd nu- 
clei, all the collective wave functions distinctly peak at positive values of (3 
deformation for all values of angular momentum. In addition, with increasing 
values of the mass number these collective wave functions get narrower, and 
their minima are located at larger values of (3. If triaxial shapes were allowed 
in the configuration mixing calculations we will expect larger mixing in the 
nuclei with softer surfaces. 

One of the best signatures of a shape transition is the behavior of the ra- 
tio between the energies of the first 4 + and 2 + states (-R4/2 = E{A + ) / ' E{2 + )) 
along the isotopic chain. This ratio varies from the value which correspond to 
vibrations around a spherical shape (-R4/2 = 2) to the characteristic value for 
excitations of a well-deformed rotor (-R4/2 = 3.33). Fig. 6(a) shows that these 
limits are experimentally fulfilled in the 144 Nd and 152_156 Nd isotopes, respec- 
tively, and also that there is a smooth transition between them. Furthermore, 
we observe a good agreement between theory and experiment, specially in 
the vibrational and rotational limits. However, we also notice that the tran- 
sition is more abrupt in the theoretical predictions than in the experimental 
results and only the 148 Nd isotope could be considered as a transitional nu- 
cleus in the calculations. Similar conclusions are extracted from the behavior 
of the reduced quadrupole transition probabilities (B(E2,2 + — > + )) shown 
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Fig. 7. Normalized excitation energies E( J^~) / E(2f) for the 144_154 Nd isotopes. 
Dots and circles correspond to experimental and theoretical values, respectively. In 
addition, U(5), X(5) and SU(3) predictions (boxes) are shown in comparison to the 
data for 144 Nd, 150 Nd and 154 Nd nuclei. Notice the different scale in the upper and 
lower panels. 

in Fig. 6(b). Experimentally and theoretically an increase of the B(E2) y s with 
increasing neutron number is observed as we are moving further away from 
the N = 82 shell closure. Again, the theoretical results display a sharper jump 
from the weak to the strong collective behavior with 148 Nd in an intermediate 
position. 

We also compare our theoretical results with the available experimental data 
(25; 15) for higher excited states in Fig. 7. The experimental energy spec- 
tra in the vibrator ( 144 Nd) and rotor ( 152 ~ 154 Nd) limits are nicely described 
with the present calculations. Additionally, we plot in Fig. 7 the predictions 
given by U(5) and SU(3) models for those cases and an excellent matching 
between both the symmetries and the experimental and theoretical results 
are observed. In the first case, the law E(J) ~ J is obeyed while in the latter 
E(J) ~ J(J + 1) is fulfilled. However, deviations from the experimental trends 
are found in 146_150 Nd. It is precisely in these nuclei where the triaxial effects 

- not included in these calculations — are expected to be more important 
(see Fig. 3). 

In Fig. 7 the predictions of the X(5) symmetry (13) are also shown. There, 
a good agreement with the experimental data for the 150 Nd isotope is found 
and, from this kind of comparisons, it has been inferred that this nucleus cor- 
responds to the experimental critical point (15). We have already mentioned 
that in our description it is the nucleus 148 Nd that looks more transitional. 
As a matter of fact, if we compare the X(5) results with the theoretical pre- 
dictions for 148 Nd, see Fig. 8, we find striking similarities in energies and 
transition probabilities suggesting that the theoretical results for the nucleus 
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Fig. 8. Normalized excitation energies and transition probabilities for the 148 Nd 
isotope calculated with the GCM method and the X(5) predictions. 

148 Nd display the X(5) symmetry. As mentioned at the beginning there have 
been also configuration mixing calculations performed in the framework of 
the relativistic model by Niksic and collaborators (20). These authors obtain 
a good agreement with the experimental results for the nucleus 150 Nd and we 
can affirm that the relativistic models also reproduce the X(5) predictions. At 
this point we are confronted with a kind of paradox, the must condition of 
Iachello to obtain the X(5) symmetry spectrum is to impose that in the crit- 
ical point nuclei the spherical and the deformed minima must coexist and be 
degenerated in order to obtain a first order shape phase transition. However, 
our calculations and the relativistic ones are able to produce X(5) results with 
potentials very much different from the one of Iachello, see Fig. 4 of the present 
work and Fig. 1 of Ref. (20). It is important to realize that the potentials used 
in the relativistic calculations as well as in ours do not satisfy the mentioned 
must condition to identify a critical point nuclei. Hence, it seems that spectra 
and transition probabilities of the X(5) type are not sufficient for the identifi- 
cation of a critical point nucleus, at least not for the type of phase transition 
described in Ref. (13). One could think that the relativistic potential of Ref. 
(20) and ours look similar and that it is by chance that the calculations re- 
produce the X(5) results. However, this is not the case because Niksic et al. 
only consider the prolate part of the potential in their configuration mixing 
calculations and neglect the mixing with oblate configurations. This is not the 
case in our calculations and as a matter of fact we have performed calculations 
both with and without the oblate part and we find large differences. As one 
can see in Fig. 9, where for the nuclei 146 Nd and 148 Nd we plot the spectra with 
and without consideration of the oblate part of the potential. In particular for 
the nucleus 146 Nd we find that in the complete calculations we are close to 
the vibrational limit and in the prolate restricted ones close to the rotational 
limit. A relevant aspect of our calculations is that they reveal the importance 
of triaxial shapes in the transitional nuclei 146 - 148 - 150 Nd and the potential 
role that they should play if considered in the calculations. One has further to 
consider that triaxial angular momentum projection of the PES of Fig. 3 will 
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Fig. 9. Normalized excitation energies E(J^) / E(2f) for (a) 146 Nd and (b) 148 Nd iso- 
topes calculated with the GCM method using only prolate (diamonds) and prolate 
and oblate configurations (circles) with the Gogny force. Furthermore, the experi- 
mental data (dots) and the X(5) predictions (boxes) are also represented. 

make them shallower, see Ref. (26), and the subsequent configuration mixing 
of triaxial shapes even more relevant. Another point which in our opinion has 
not received enough attention in earlier works is the presence of negative par- 
ity states with J 71 = 3~ and J w = 1~ experimentally observed in the low-lying 
energy spectra for the 146_150 Nd (25; 27) and Samarium isotopes (25). These 
data indicate the relevance of the octupole degree of freedom in this region. 
Though our present parity conserving calculations ignore this degree of free- 
dom, earlier calculations with the Gogny force which had considered it, have 
shown (28) that the Nd isotopes are rather soft against octupole deformation. 
Furthermore, these calculations were able to reproduce these negative parity 
states (28) up to high angular momentum. 



In conclusion we have studied the spherical to prolate deformed shape transi- 
tion in the Neodymium isotopes ( 144_154 Nd) with state-of-the-art beyond mean 
field calculations with the Gogny force. We have applied the generator coor- 
dinate method with particle number and angular momentum projected axial 
wave functions to this problem. We have found a good agreement with the 
experimental data, specially in the vibrational and rotational limits, although 
the calculations are not able to reproduce quantitatively the transitional re- 
gion. As a matter of fact, the transition is theoretically more abrupt than in 
the experiments. Concerning the transitional region we would like to underline 
the following points: First, there are several intrinsic potentials which display 
spectra and transition probabilities corresponding to the X(5) symmetry and 
which do not satisfy the condition, given in Ref. (13), to identify them with 
the critical point of a first order shape phase transition. Second, the simplistic 
view point of considering in a phase transition all degrees of freedom frozen 
but the axial quadrupole moment is not supported by microscopic calculations. 
On the contrary, these calculations seem to support a general relaxation of the 
most relevant degrees of freedom. Third and last, the consideration of the beta 
deformation as the order parameter of a hypothetical shape phase transition is 
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not well founded. All these points question the interpretation of nuclear shape 
changes as nuclear shape phase transitions. The final interesting task for the 
future is to see whether the inclusion of additional degrees of freedom (spe- 
cially the triaxial one) leads us to a quantitative description of the transition 
region. 

This work has been supported in part by the DGI, Ministerio de Ciencia y 
Tecnologia, Spain, under Project FIS2004-06697. 
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